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RESONANCE TYPE INSTABILITIES 
IN THE GASEOUS DISKS 
OF THE FLAT GALAXIES 
HI. The gyroscopical resonance type instability 

C.M.Bezborodov, J.V.Mustsevaya and V.V.Mustsevoy 

The stability of thin homogeneous gaseous disk with step discontinuity of velocity and den¬ 
sity profile has been investigated analytically in the limitary case of large compressibility and 
numerically at spatially separated velocity and density kinks. The expressions for reflection and 
transmission coefficients at the axisymmetrical vortex sheet are derived. The possibility of gyro¬ 
scopical resonance type instability development in addition to Kelvin-Helmholtz and centrifugal 
instabilities is shown and its physical mechanism is clarified. It has been discussed the possibility 
of application of obtained results to the stability analysis of the gaseous disks of the real flat 
galaxies. 

Introduction. 

The applied problems of astrophysics about stability of supersonic currents (i.e. disks 
of accreting matter, galactic disks, jets from young stars etc.) from the academic point of 
view are interesting hrst of all by allowing to allocate many specihc effects and mechanisms 
of amplihcation of disturbances in the pure state. This circumstances are related to absence 
of such complicating factors as heat exchange with hrm borders, non-absolute rigidity of 
these borders, and also extremely small viscosity. Therefore it is possible to combine 
successfully relative simplicity of models with practical application of results obtained 
at their analys22-08-98 09:39pm22-08-98 09:39pmis. In the present part of work with 
reference to gaseous disks of flat (spiral) galaxies amplihcation of unstable disturbances in 
a waveguide layer between jumps of rotation velocity and density caused by over-rehection 
is discussed. 

The effect of approaching inhnity of rehection and transmission coefficients of a hat 
monochromatic sound wave falling at certain angles on a hat surface of supersonic vortex 
sheet (so called the ehect of over-rehection) was discovered by Miles (1957) and Ribner 
(1957). Subsequently the similar ehect was investigated for magnetohydrodynamic jump 
(Acheson 1967). It was also shown that the presence of jump is not obligatory for over- 
rehection and that it also takes place on a critical layer in a how with smooth velocity 
distribution (Kolyhalov 1985) (that is on a layer where how speed and wave phase speed 
along how speed coincide). The mechanism of over-rehection is caused by that the energy 
hux in a transmitted wave is directed to jump (Miles 1957, Ribner 1957). Thus the 
transmitted wave outgoing from jump (or critical layer) picks up kinetic energy from the 
basic current and transfers it to a rehected wave through jump (see e.g. Landau & Lifshits 
1988). 

Although Betchov & Criminate (1967) marked that if to arrange a rehecting surface 
(acoustic screen) in parallel to a critical layer, an acoustical resonance type instability may 
develop. However, their speculations did not go further than statement of such opportunity. 
The hrst rather serious researches of acoustic resonance type instability are probably works 
of Ferrari et al. 1982, Payne & Cohn 1985, Hardee & Norman 1988 (see also Norman et al. 
1982), in which stability of astrophysical jets in models of a hat layer and cylindrical jet 
with discontinuous borders was investigated. The large number of works on the specihed 
theme have appeared subsequently but all of them concerned to plainparallel gas hows. 
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The unstable modes caused by over-reflection from a vicinity of corotation in gaseous 
disks with monotonous dependence of rotation velocity on radius were found for the first 
time by Papaloisou and Pringle (1985, 1987). Then such modes (and also other resonant 
modes, for which amplification a continuity of rotation velocity is essential) were inves¬ 
tigated by Glatzel (1987), Savonije & Heemskerk (1990), Papaloizou & Savonije (1991) 
etc. Nevertheless as for disks with continuous radial parameters distribution the linearized 
equations of gas dynamics can be investigated only numerically, there was the number of 
unsolved questions connected to dependence of characteristic time of development of un¬ 
stable disturbances and angular phase rotation velocity of a pattern created by them upon 
unperturbed disk parameters. Besides, in quoted works the consideration was performed 
with reference to accretion disks. 

In the present part of the article the effect of over-reflection and instability caused by it 
in the elementary discontinuous model of a thin gaseous disk is considered. Its applicability 
is caused by observable peculiarities of density and rotation velocity distributions in disks 
of flat galaxies. 

The basis for consideration is such fact that the wave can be many times reflected from 
rather sharp density jump and from corotation radius, and from the latter with strengthen¬ 
ing (Savonije & Heemskerk 1990). It is clear that even when “centres” of smoothed jumps 
of velocity and density distributions coincide, specified points of turn generally speaking 
are separated on radius. Thus in a disk the waveguide layer presents, in which wave collects 
energy in due course. 

Is the most simple model of such waveguide in a gaseous disk is model of two jumps, 
namely vortex sheet^ and internal rather it contact jump. The choice of discontinuous 
model for the analysis requires some explanations, as to the present there is the number of 
works where stability of gaseous disks with continuous velocity and density profiles (see, 
for example, Papaloizou & Pringle 1985, 1987, Glatzel 1987, Savonije & Heemskerk 1990, 
Papaloizou & Savonije 1991, Torgashin 1986, Morozov 1989, Morozov et al. 1992) was 
analyzed. 

However in a disk with continuous parameters distributions the instability as a rule 
is developed simultaneously by several mechanisms (for example, over-reflection and over¬ 
transmission of waves in a corotation vicinity (Savonije & Heemskerk 1990), resonant 
radiation of waves from corotation radius (ibidem), centrifugal (Torgashin 1986, Morozov 
1989, Morozov et al. 1992), centrifugal-resonant (Morozov et al. 1992), Kelvin-Helmholtz 
instability (KHI) mechanism (Torgashin 1986) etc.). So it is almost impossible to tell them 
from each other and to reveal their relative influence. In the contrary the discontinuous 
model will allow us in section 3 to show that instability developing in it is of a really 
new type. Though the physical mechanism of this instability has similar features with 
the mechanism of acoustic resonance type instability in plainparallel current (Payne & 
Gohn 1985), there is conceptual difference consisting that a resonance occurs rather at 
sound frequency than at gyroscopic one. Accordingly frequency areas of unstable modes 
localization are determined not by “sound” frequency |k|cs but gyroscopic one f2in;ea;(?^± 
2), where m is azimuthal mode number (number of an arms of spiral pattern). Besides 
such model will allow to allocate in the pure state unstable superficial modes: centrifugal 

^ In discontinous model corotation radius with necessity coincides with vortex sheet, because 
phase speed of a unstable wave should lay between the minimum and maximum speeds of the 
basic current, that it is possible to prove with the help of a theorem (Kolyhalov 1985). 
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connected with velocity jnmp and mode connected to density drop, and to investigate 
dependence of their properties on parameters of a disk. Thns an independent problem 
is to determine the expressions for wave falling on axisymmetic vortex sheet reflection 
and transmission coefficients (section 2), as thongh they have large importance for onr 
consideration bnt were not received former. 

Finally, it seems nsefnl to discnss in section 4 the probable astrophysical applications 
of obtained resnlts. 

In the present part of the work we shall nse the references to the formnlae from Parts 
I and II, adding in this case a symbol “1.” or “2.” accordingly before nnmber of the 
formnla. 

1. Model and dispersion equation. 

The basic characteristic featnres of model nsed in the present part of onr work are the 
same as in model described in a Part II. We will therefore stop only at differences. 

The model contains not one bnt two jnmps: contact jnmp at radins Rp and vortex 
sheet at Let’s assnme for dehniteness Rp < Rq as it is seems snch sitnation is more 
characteristic for real objects (at least it is valid for the Galaxy). The radial dependences 
of nnpertnrbed parameters are set by relations: 


Po(G z) = Pin{z) + iPex{z) “ Pin{z)) 9{r - Rp), 
Cs{r,z) ^ Csir,{z) |1 + 


Pin{z) 


Pex{z) 


l]0{r- Rp) 


0(r, z) = Qin{z) + {Qex{z) - ^in{z)) 9{r - Rq), 
Poir, z) = c^(r, z)po{r, z)/'f ^ const{z). 


( 1 ) 


As the magnitndes of distnrbances shonld be limited in zero and at inhnity, the solntion 
of the appropriate modihed Bessel eqnation in the given model takes a form: 


BImi,kme^) R K.uiikme^) i ^ ^ Rq^ 

DKm{k ), ^ > Rn- 

Here kin and kex are dehned by (2.6) and in area between jnmps 

- mftin)^ 

^me 2 

ex 

It is obvions the conditions at jnmps do not differ from (2.8), (2.9): 

p{Rp + 0) - p{Rp - 0) =0, 
p{Rfi + 0) - p{Rq - 0) =Pex^iRn) Rfii^in - ^Ix): 

^{Rp + 0) - ^{Rp - 0) =^(i?n + 0) - ^{Rq - 0) = 0. 


( 2 ) 


( 3 ) 


( 4 ) 


Snbstitnting the solntion as (2) in conditions (4) we come to the system of linear algebraic 
eqnations on coefficients A, B, C, D. A condition of simnltaneity of this system will be 
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equality to zero of its determinant. Excepting neutral gyroscopic modes of fluctuations 
u) = (m ± 2)Qin, u) — (m ± 2)flex we finally write: 


1 

1 

1 

0 



Qai^l 


0 

- a^ex [1 

0 

6k 

6i 

Aq^ — {x — mqfl - 

0 

dROime 


[4 - {x 

-m)2)] 


The equality (5) represents the dispersion equation for disturbances in considered model. 
Designations of Parts I and II in combinations and setting by (1.8) and (1.9) are 
here kept; top index specihes that they are written on radius of velocity or density jump, 
the bottom index takes values “in”, “ex” or “me” (with the account = f^m)- Besides 
following designations are introduced: 


5r = 


.fm (^me Ra), 
ImikmeRp) 


5k = 


Km{,k me Rn) 

Kmi,k me Rp) 


( 6 ) 


Before solving directly the dispersion equation (5) in the following section we con¬ 
sider a problem on reflection of a neutral wave from axisymmetic vortex sheet to have an 
opportunity to apply its results at interpretation of the instability mechanism. 


2. Reflection and transmission coefficients. 

From theoretical works (Miles 1957, Ribner 1957) it is known that at incidence on 
plainparallel vortex sheet the neutral wave can reflect and refract with amplihcation on 
amplitude (over-reflection effect). At presence of the second reflecting surface (for example, 
density jump or kink) spatial resonator will be formed, in which the wave can collect energy 
in due course. Therefore the over-reflection is capable to result in development of insta¬ 
bilities of a waveguide-resonant type, in particular gyroscopic resonance type instability 
considered in the given section. 

For plainparallel shift currents the problem of over-reflection up to now is investigated 
well even for systems with magnetic helds (Acheson 1967). However, coefficients of reflec¬ 
tion and transmission for a neutral wave falling on axisymmetic jump were not obtained 
till now. 

Note at once that in cylindrical geometry there is the allocated direction on centre. 
It leads to distinction of coefficients for waves falling on axisymmetic jump from within 
and from outside. Nevertheless, in a limit of large jumps radii or azimuthal short-wave 
disturbances the expressions for coefficients received in our model should lead to similar 
expressions from a plainparallel problem. 

Following Miles (1957) and Ribner (1957), we shall consider incidence of a neutral 
monochromatic wave on vortex sheet being not interested in this section by a question of 
its stability (on formal correctness of such approach see Landau & Lifshits 1988). 

Further in this section the index “i” concerns to a incident wave, “r” — to reflected, 
index “t” — to transmitted ones. The magnitudes in internal area from jump are desig¬ 
nated by index 1, outside — 2. 

a) Incidence from within on density and velocity jump 
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By applying results of a Part II and condition of limitary absorption (i.e. the distur¬ 
bances should fade as they propagate) at once write out peak functions of pressure of all 
three waves: 

Pi = AKra{kir)] Pr = S/^(/cir); pt = CKm{k 2 r). (7) 

Thus displacement magnitudes on dividing border are 


6 


— 


6 = 


AKmikiR)ai -|- BIm{kiR)l3i 

pic^^kiR 

CKm{k2R)o'2 

P2C^2^'^R 


( 8 ) 

(9) 


The parameters cti, (3i are governed by the formulae (1.8), (1.9). 

Now introduce complex reflection and transmission coefficients with respect to pres¬ 
sure by a usual way: 


_ Rlmi^kxR^ ^ _ Pt CKui(^k2R) 

Pi AK^{kiR)' Pi AKmikiR)' 


( 10 ) 


Taking into account general conditions (2.8), (2.9) come to the equations on 


klR^ a2ji-) _ A ^(-) _ . 
klR‘^ ai ai 





M2(l + g) a2 

2p.‘^ klR'^ 






Solving this system let us write coefficients in an obvious kind: 


(11) 

( 12 ) 


^(_) Qk\R^a2 — p^k'^R'^cti + |(1 + ( 5)(1 — q^)M‘^aia2 
^ ^ p^kjR^Pi - QklR^a2 - i(l + Q)(l - q^)M^/3ia2 ’ 

ji-) ^_ Qk^R^iPi - ai) _ 

iPklmPi - QklRPai - 1(1 + Q)(l - q^)M^Pia-2' 


(13) 

(14) 


b) Incidence from outside on density and velocity jump 

In a case when the neutral wave falls on vortex sheet from the large radii the solution 
of modihed Bessel equation takes a form: 


Pi = AImik 2 r)- Pr = BK^{k 2 ry, pt = CI^{kir). 

With the help (15) write down displacement on jump: 

CIm{kiRp)(R 


6 = 6 
6 = 6 + 6 


PicliklR 

AIm{k2R)/32 + BKm{k2R)c^2 
p2cl2klR 


(15) 


(16) 

(17) 
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Applying again the condition on jnmp and making simple transformations we obtain co¬ 
efficients and 


- QklR^P2 - i(l + Q)(l - q^)M^l3il32 
QklR^ai - ^‘^k^R^Pi + i(l + Q)(l - q^)M^l3iai ’ 

j{+) ^_ ^‘^kfR^jai - (32) _ 

QklR^^ai - fRklR^^i + i(l + Q)(l - q^)M^(3iai' 

Note that classical connection between TZ and T (7= 1 -f7^) known for a plainparallel 
problem (Miles 1957, Ribner 1957, Landan & Lifshits 1988) in onr case is valid on density 
jnmp only as jnst on it the bonndary conditions coincide identically with those of a fiat 
problem, namely both ^ and p are continnons. 


(18) 

(19) 


It is interesting to note taking the determinant (5) it is possible to write dispersion 
eqnation throngh the reflection coefficients prodnction: 


n 


(-|-)^( —) Im{.kmeRvi) ^mikmeRp) 

ikme^p) K^(^k me Rn) 


( 20 ) 


Here the index “p” marks reflection coefficient on density jnmp {R = Rp, q = 1) and 
index “fl” — on velocity jnmp {R = Rq, Q = p = 1). The latter eqnation allows 
to demonstrate evidently that reflection coefficients are associated with the instability 
mechanism that we consider in the following section. Really, in the limit of short radial 
wavelength {\kmeRp\ ^ w) from (20) we find: 


^ - exp {2kme(Rn - Rp)} ^ exp 


-2iD — {x 

p 



Dimensionless width of a backlash between jnmps is here intodnced: 

jj _ R^ - Rp 

Rfi 

Taking a sqnare of the (21) modnle we get: 




~ exp IaD—I mx \ = exp 




■'Sex 


( 21 ) 


( 22 ) 


Thns for nnstable {Imu > 0) distnrbances 
follows 




(-) 

n 


> 1 is always valid. 


From (22) 


Imu ~ 


—In 

2Tr 


K<+'K<i > 


(23) 


where Tc is characteristic rnn-time of a sonnd wave between jnmps. Note that the resnlt 
(23) is typical for models admitting development of aconstic resonance type instabilities 
(Payne & Cohn 1985). 


3. The stability analysis. 
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Before receiveing the approximated solutions of the equation (5) it is necessary to make 
one remark. For existance of over-reflection and unstable reflective harmonics caused by 
it a nonzero wave energy flux in a radial direction is necessary at least in two areas (at 

r > Rq and at Rp < r < Rq) the same as and in a flat case (see Morozov et al. 1991). It 

is possible if the eigenfunctions p and ^ in these areas are oscillating on radial coordinate 
(as for value of group velocity a value of radial component of wave vector kr is crucial). 
The latter may occur only when arguments of modihed Bessel functions are imaginary (the 
preceeding concernes to a limit Imu) 0, as for such disturbances only an energy flux is 
dehned). Writing k‘^ = —k‘^ in each area we come to restriction on a real part of frequency, 
namely, the following inequalities should not occur: 

(m — 2)Qin < Reu) < {m + 2)Qin, (24) 

(m — 2)flex < Reu) < {m + 2)flex- (25) 

As for increasing disturbances 

mflex < Re oj < mVtin ( 26 ) 

is valid (at Vt^x < the reflective harmonics can take place only in following frequency 
range: 

(m -|- 2)flex < Reuj < (m — 2)0^^, (27) 

that it is possible at m > 2. Note that here is the strict inequality, as neutral gyroscopic 
modes have already been excluded from consideration. 

Will jump-orthogonal component of a wave vector be real or imaginary (and accord¬ 
ingly presence or absence of energy flux in this direction) is determined in a flat case by 
a relation of squares of Doppler w — k||V and sound k\\Cs frequencies, and in considered 
case — by Doppler u — mVti and gyroscopic 2Vti ones. And in the latter case changing Cg 
simply scales /c (|/c| —> cx) at Cg —0 and |/c| —0 at Cg —cx)), not breaking proportions 
between Imk and Rek. Accordingly the resonance will take place on harmonics of gyro¬ 
scopic frequency (m — 2)Vtin in a layer of gas between velocity and density jumps. This 
allows on the one hand to limit ourselves in searching of reflective harmonics by permitted 
area of frequencies (27) (this area is limited by continuous lines on Fig. 3.1) and on the 
other hand to set the value (m — 2)0^^ as initial approximation of frequency of the basic 
unstable mode in a waveguide layer between jumps. 

Let us obtain the solution of the dispersion equation (5) in a limit of strong compress¬ 
ibility of media in all three areas divided by jumps (M = RciQin/csin 3> 1, M/p ^ 1). 
Remember, magnitude of arguments of modihed Bessel functions is determined in the ba¬ 
sic by Mach number. It allows to put \kiRp\ 3> m; in this case the dispersion equation is 
simple to bring to a kind: 

Qilva + p,^ x-m 

ptha + Q —{x — mq)+i^{l — q‘^)^ 

where a = kmeRfi D = a/ 4 — (a: — m)^. For the case of uniform on density at ^ = const 

disk (Q = p = 1) from (28) result (2.12) describing centrifugal instability mode directly 
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follows. In general case {Q 7 ^ 1) if only D ^ 0 has not place with accnracy to the terms 
of order — //)/(Q + /u) | 1 we hnd: 

x ~-m(l + q) + i —(1 - (29) 

^ ^ Uj 

Expression (29) means that if the jnmps are not too close the centrifngal mode of velocity 
jnmp in the main order is not “sensitive” to the density jnmp. 

For the fnrther calculations transform (28) as follows: 


Q{x — m) + fx{x — mq) — iM{l — q^) 
n{x — m) + Q{x — mq) — — q^) 


(30) 


In a case of very close incoincident jumps {D 1) the left-hand part of (30) can be 
presented as: tha ~ —iMD{x — m)/fx. Substituting this result in (30) we get the square 
equation on x, but despite its simplicity the solution appears extremely inconvenient for 
the analysis because of bulkness. Therefore let us limit by a case D ^ 0 ai which we hnd: 

x ~— ^-—{m{Q + qfj.)+iM{l-q'^)}. (31) 

H + Q 

It is easy to note, that (31) does not coincide with the asymptotic solution (2.11) for 
concurrent jumps. Thus the limiting transition to a case T> = 0 is absent. It is a direct 
consequence that we work with model of “jump” with hnite thickness (^ An, Ap, where 
An, Ap are characteristic radial scales of jumps of angular velocity and density), in which 
as it was specihed by Fridman & Khoruzhii (1993) the large role is played with a specihc 
kind of equilibrum parameters inside “jump” prohle. It is obvious, the result (31) could 
be obtained also in model of concurrent jumps if to set “asymmetrical” concerning middle 
of “jump” p(r) and Q{r) distributions that would change boundary conditions (Fridman 
& Khoruzhii 1993). The physical sense of (31) becomes perfectly transparent if to write 
growth rate (31) in the dimensional form: 


_ PexRaj^in 

CsinPin T CsexPex 

So the converce characteristic time of instability development is determined by the ratio of 
centrifugal force density drop on angular velocity jump to the sum of media wave resistances 
on both sides from a layer containing both jumps. 

To hnd out rehective harmonics we transform the left-hand part of (30) by presenting 
dimensionless frequency as x = Xq + Sx and |ha:| Xq: 



th 


^ MD ~itg 


MD 

P 


a/ (xo — m)2 — 4 + C 5x 


(33) 


MD (xo — m) 

P a/(xo - m)2 - 4 
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where 



Transforming (33) fnrther, we assnme a/(^cq — mY — 4 = 7 rn. 

Then (33) degenerates to ttg( 7 rn + C' 5x) = ttg(C' 5x). Making the additional assnmp- 
tion \C 5x\ -C 1 we get tg(C' hx) ~ C 5x^ with respect of told let ns write approximated 
expression for freqnency of reflective harmonics (nnmber of harmonic n = 1, 2, 3,... deter¬ 
mines nnmber of zeros of eigenfnnctions between jnmps): 


xq = m 


4+r^v- 

\md) ’ 


dx ~ 


m{Q iiq) - xo{Q + n) + iM{l - q^) 

(Q + /u) + fi^(l - q^) -I- iC [xo{Q + n) - m{n + qQ)] 


Remind, by virtne of told in the section beginning, xq lays in limits: 


(34) 

(35) 


(m -I- 2)q < Xq < (m — 2). (36) 

At last let ns hnd approximated expression for basic (with n = 0) nnstable wavegnide 
mode. Note, by simple snbstitntion n = 0 in (34), (35) the freqnency of this mode cannot 
be described as then C tnrns to inhnity. 

We believe x ^ m — 2. Then in approximation of \kinRp\^ \kmeRQ.\ -C 1, \kexRQ.\ 3> 1 
snbstitnting freqnency as x — m — 2 + 5x in the eqnation (5) we determine: 

^_ 3[l-g" + S(m- 2 -m 9 )] _ 

ill-E]-^ + |(l+e)+=](l-42 + ]S(m- 2 -m,))’ 

where S = (1 - g)(l - Df^. 

Remarkable pecnliarity of expressions (34), (35), (37) is the fact that they describe 
nnstable roots both at > Vtex and at Vtex > klin and in an essentially snpersonic case 
{Ru\^in — ^ex\/Csex ^ 1 )? and are similar to analogons asymptotics of reflective harmonics 
and basic mode developing between density and velocity jnmps in a plainparallel flow (see 
Morozov et al. 1991). According to told in the beginning of section instability described 
by expressions (34), (35), (37) we shall name the gyroscopic resonance type instability 
(GRTI). The absence of centrifngal snppession of instability at Vtex > ^in is cansed by 
over reflective character of its development mechanism that is by specihcally wave effect 
not connected directly to mass forces density distribntion in a disk. 

The case of small compressibility (M 1) cannot be investigated analytically becanse 
of bnlkness of the approximated algebraic eqnation obtained in this limit. 

The resnlts of the nnmerical solntion of (5) (see Fig. 3.2-3.7) as a whole conhrm 
conclnsions made on the basis of asymptotic analytical research and at essential compress¬ 
ibility (Rn(f^m — ^ex)/csex ^ 1) ^md Small relative distance between jnmps D <^1 these 
resnlts show presence besides mode of centrifngal instability a discrete set of poorly nn¬ 
stable GRTI modes. The phase angnlar velocity of rotation of a pattern of these modes 
Vtp = Reu/m extremely weakly depends on Q and q (Fig. 3.2, 3.4) and essentially changes 
with D and M change (Fig. 3.5, 3.6), remaining thns in permitted freqnency area (27) (see 
also Fig. 3.1). Exception makes basic (not having zeros of eigenfnnctions between jnmps) 
GRTI mode for which in all parameters ranges Reu nnsignihcantly exceeds (m — 2)Qin- 
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The number of reflective (n > 1) harmonics is increased with growth of Mach number, 
azimuthal mode number m and distance between jumps reduction. 

As it is seen from Fig. 3.8 the asymptotic solution of (34), (35) gives only qualitative 
conformity to numerical results. It is connected that at approximated taking square roots 
in kiRci we have neglected a square of gyroscopic frequency in comparison with a squared 
Doppler one, whereas the effect of over-reflection resulting in GRTI occurrence is caused 
by a resonance just between these frequencies. Nevertheless, (34), (35) begin to work well 
at large numbers m when the calculation of modihed Bessel functions on the computer is 
complicated. 

Let us specify at last, that as well as in a Part II the considered disturbances with 
Vz = 0 are allocated by a natural way by the instability mechanism itself. Really, in this 
sense all told in Part II is valid for centrifugal instability mode, and resonant harmonics 
with a wave vector componenet kz comparable with kr and k^p during the time of h/cg 
order (here h is the characteristic disk half-thickness) leave a waveguide layer. As this time 
is near to inverted growth rate of resonant modes, such disturbances will have no time to 
amplify up to appreciable amplitudes. If to take into account an opportunity of returning 
of these disturbances to a waveguide layer because of refraction, their growth rate anyway 
will be less than at disturbances not having kz- 

4- Conclusions. 

Let us generalize numerical and analytical results by formulating the conclusions about 
stability of a uniform on pressure at ^ = const gas disk containing rotation angular velocity 
jump and internal rather it non-homentopic density jump (in models of “jumps” of hnite 
thickness, that is ^ A^, Ap, where Aq, Ap are characteristic radial scales of angular 

velocity and density jumps). 

1. In a limit of weak compressibility in considered system only one unstable mode 
resulting Kelvin-Helmholtz instability development on vortex sheet can develop. 

2. With growth of compressibility the development mechanism of the vortex sheet 
superhcial mode essentially changes from Bernoulli effect resulting in Kelvin-Helmholtz 
instabilities up to centrifugal one in a supersonic case. 

3. The presence of density jump in the system appreciably impacts to centrifugal 
instability development only when the jumps are very close, i.e. D 1. 

4. At essentially supersonic velocity difference on jump and m > 2 a discrete set of 
weakly unstable modes appears. These modes differ by number of zeros of eigenfunctions 
between jumps. The growth rate of these modes appreciably differs from zero only when 
D 1 and their amount grows with growth of compressibility on jump, azimuthal mode 
number, with distance between jumps reduction and practically does not depend on other 
parameters. These modes are caused by over-reflection on velocity jump and collect energy 
in a waveguide layer between jumps exponentially in time, picking it up from kinetic energy 
of the main movement (rotation) of outer gas layer. 

5. The effect of over-reflection is caused by a resonance between Doppler and gyro¬ 
scopic frequencies^ and the waveguide unstable modes described above are displayed near 
to the basic gyroscopic frequency and its higher harmonics. So the given instability is 
named as gyroscopic resonance type instability (GRTI). 

^ The last in considered model with discontinuous velocity profile represents a degenerate case 
of epicyclic frequency. 
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6. As against centrifugal instability the GRTI develops and at “negative” velocity 
drop on jump (when the periphery of a disk rotates faster than the central region). 

Finishing the article, we make some hnal remarks on probable applied meaning of 
obtained results. 

It is necessary to note, GRTI by virtue of small scale and small growth rates is hardly 
capable to compete to centrifugal and Kelvin-Helmholtz instabilities as the generator of 
large-scale spiral structure, its role rather can be reduced unless to creation of multi-arm 
secondary (imposed on basic) pattern (see Morozov et al. 1992). At the same time it is 
known, the similar resonant Papaloisou-Pringle instabilities and acoustic resonance type 
ones are much more weakly subject to stabilizing influence of velocity jump smoothing 
and besides are multimode (Savonije & Heemskerk 1990, Papaloizou & Savonije 1991, 
Mustsevoy & Hoperskov 1991). Add that the resonant modes by virtue of its development 
slowness do not result to rough rearrangement of system parameters distributions and 
to system exiting on stability margin. Therefore the simultaneous development of large 
number of various harmonics of different azimuthal GRTI modes probably is capable to 
lead to effective removal of an angular moment from a vicinity of a waveguide layer on disk 
periphery and to smooth durable change of radial profiles of density and rotation velocity. 

Make at last one more remark on possible GRTI and other hydrodynamic resonant 
instabilities influence on magnitude of stars velocities dispersion in a galactic disk. 

Besides gyroscopic resonance type instability considered in the previous sections the 
development of other resonant instabilities (e.g. Papaloisou-Pringle instability (Papaloizou 
1985, Savonije & Heemskerk 1990), resonant-centrifugal one (Morozov et al. 1992)) is 
possible in a differentially rotating gaseous disk. All of them are multimode and are 
characterized by discrete sets of radial, azimuthal and transverse to a disk wave numbers 
and frequency eigenvalues. Thus in result of development of resonant type instabilities 
in a gaseous galactic disk there is the complex system of disturbances of velocity and 
thermodynamic parameters with hierarchy of spatial and temporary scales. It is rather 
essential that the energy of these disturbances exceeds sound energy or tends to it from 
above. The basis for the last statement give laboratory (Norman et al. 1982) and nonlinear 
numerical (Norman & Hardee 1988) modelling of supersonic jets. In these experiments it is 
shown that the growth of disturbances caused by resonant instabilities stops at a nonlinear 
stage at the formation of system of weak slanted shock waves (disturbance with wavelength 
comparable to characteristic jet size) on a background of advanced acoustic turbulence (the 
most short wavelength disturbances). On the other hand it is known that in galactic disks 
rich by gas (~ 10%), the gravitational stability of a disk is determined by more “cold” 
gas componenet (Jog & Solomon 1984). Therefore in numerical experiments on modeling 
binary (gas -|- stars) disk the gravitational instability of a gas subsystem develops and 
results in formation of giant molecular complexes (GMG). In turn relaxational processes 
(stars scattering on massive GMG) result to fast “warm up” of a star disk. As a result stars 
velocity dispersion becomes so high that the star disk appears gravitationally stable even at 
a unstable initial condition (described scenario was observed, for example, in experiments 
of Shlosman & Noguchi 1993). So in numerical experiments the specified process could 
result in a conclusion of impossibility of spiral structure or bar generation by development 
of gravitational instability of a star disk. 

As the margin of gravitational stability is determined by balance of energy of an 
ordered movement under action of gravitational forces and stochastic thermal movements 
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and as the energy of small-scale system of weak slanted shock waves cansed and snpported 
by mnltimode resonant instabilities, dissipates in thermal becanse of viscosity (thongh 
small), it is necessary to make a conclnsion that the really gas disk in an extended vicinity 
of a wavegnide layer is warmed np and conseqnently appears to be more gravitationally 
stable. Besides velocity and density flnctnations cansed by resonant instabilities shonld 
slow down GMC formation, partially to destroy already formed complexes and accordingly 
make scattering of stars on them less effective. 

So the acconnt of the really working resonant type instabilities onght canse rednction of 
stars velocities dispersion in nnmerical modeling “gas -|- stars”. We have nsed a snbjnnctive 
inclination as the resonant modes conceptnally are short-wave and cannot be observed 
in really feasible experiments as their characteristic scales are less or comparable to the 
characteristic size of a cell and they are inevitably absorbed by nnmerical viscosity. 

Acknowledgement. One of ns (VVM) is gratefnl to INTAS for snpport of this work by 
grant project N 95-0988. 
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Fig. 3.2. Dependence of dimensionless freqnency (a) and dimensionless growth rate 
(b) from relative density drop on internal jnmp. The letter designations have dispersion 
cnrves of centrifngal mode, the digits correspond to nnmber of harmonics of gyroscopic 
resonance type instability. Centrifngal instability growth rate is not shown since it signif¬ 
icantly exceeds GRTI growth rate. 
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Fig. 3.5. Dependence of dimensionless freqnency (a) and dimensionless growth rate 
(b) from Mach nnmber. 
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Fig. 3.6. Dependence of dimensionless freqnency (a) and dimensionless growth rate 
(b) from the ratio of density jnmp radins to velocity jnmp one A = Rp/R^ = 1 — 11. 
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Fig. 3.8. Comparative behaviour of dispersion curves by results of the numerical 
solution of (5) — continuous lines and on asymptotics (34), (35), (37) — dashed lines. It 



is visible concurrence only qualitative. 





